THERE ARE SALEM NUMBERS OF EVERY TRACE 



JAMES MCKEE AND CHRIS SMYTH 

Abstract. We show that there are Salem numbers of every trace. The nontrivial part 
of this result is for Salem numbers of negative trace. The proof has two main ingredients. 
The first is a novel construction, using pairs of polynomials whose zeros interlace on the 
unit circle, of polynomials of specified negative trace having one factor a Salem polynomial, 
with any other factors being cyclotomic. The second is an upper bound for the exponent of 
a maximal torsion coset of an algebraic torus in a variety defined over the rationals. This 
second result, which may be of independent interest, enables us to refine our construction 
to avoid getting cyclotomic factors, giving a Salem polynomial of any specified trace, with 
a trace-dependent bound for its degree. 

We show also how our interlacing construction can be easily adapted to produce Pisot 
polynomials, giving a simpler, and more explicit, construction for Pisot numbers of arbi- 
trary trace than previously known. 



1. Introduction 

A Salem number is an algebraic integer greater than 1 whose other conjugates all lie in 
the closed disc \z\ ^ 1, with at least one on \z\ = 1. Our main result is the following. 

Theorem 1. For every negative integer —T there is a Salem number of trace —T and 
degree at most exp exp (22 -|- 4T log T) . 

It is easy to produce Salem numbers of any nonnegative trace, so the title of the paper is 
justified. The interest in this result is that, until now, all Salem numbers found had trace 
no smaller than —1 ( JH])- Furthermore, it is now known that a Salem number of degree 
d ^ 10 has trace at least [1 — (i/9j , and it seemed conceivable that there was a finite lower 
bound for the trace. For more details see the end of the paper. 

To provide a little background, we give a brief sketch of some facts about Salem num- 
bers. The minimal polynomial P{z) of a Salem number r is reciprocal, that is, it satisfies 
z'^'^^^P{l/z) = P{z), so that is a conjugate of r, and the coefficients of P are "palin- 
dromic". All conjugates of r apart from r and lie on \z\ = 1, and P{z) has even degree. 
For every e > and Salem number r there is a A G Q(t) such that for A; = 0, 1, 2 . . . all Xt'' 
have distance at most e from an integer. If a number field K contains a Salem number r 
of full degree [K : Q] then every full degree Salem number in K is a power of the smallest 
such Salem number in K. It is not known whether there are Salem numbers arbitrar- 
ily close to 1. If "Lehmer's conjecture" is true, then there are not. The smallest known 
Salem number 1.176280818 ■■■ , discovered by Lehmer in 1933, has minimal polynomial 
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L{z) = z^^ + — z'^ — ~ — — + z + 1. The polynomial L{—z) had just appeared 
(in 1932) in Reidemeister's book as the Alexander polynomial of a pretzel knot. For 
recent connections with knot theory, see E. Hironaka P^. The polynomial L{z) can also 
be obtained from the characteristic polynomial Eio{x) of the (adjacency matrix of the) 
Coxeter graph Eiq = ' ' I by the transformation L(z) = z^Eiq{z^^'^ + z~^/'^). There 
are currently 47 known Salem numbers less than 1.3 (see Mossinghoff |llj). 

Salem numbers are closely related to Pisot numbers, which are much better understood. 
A Pisot number is an algebraic integer greater than 1 whose other conjugates all lie in the 
open disc \z\ < 1. For every Pisot number 6, the distance of 6'" from the nearest integer 
tends to as n — s> oo. The set of all Pisot numbers is a closed subset of the real line. Every 
Pisot number is a limit point of Salem numbers, and Boyd jU p. 327] has conjectured that 
Pisot numbers are the only limit points of Salem numbers. If this is true, then the set of 
all Pisot and Salem numbers is also closed. See Bertin et a/ [Ij, Boyd Ghate and 

Hironaka jTj and Salem |T7j for these and other results about Salem and Pisot numbers. 

It is already known (see ^^E]) that there are Pisot numbers of every trace. However, 
we can greatly reduce the known upper bound for the smallest degree of a Pisot number 
of given negative trace. 

Theorem 2. For every negative integer —T there is a Pisot number of trace —T and degree 
at most the sum of the first 2T + 4 primes. 

This sum is asymptotic to 2T^\ogT. The simple examples z^ — z — 1^ z"^ — z — 1 and 
z — n{n ^ 2) of minimal polynomials of Pisot numbers then show that there are Pisot 
numbers of every trace. 

Computations for negative trace down to —25 (see Section ^ indicate that the upper 
bound on the degree in Theorem ^should be comparable with that in Theorem|21 However, 
a proof of this does not seem within reach at present. 

In for infinitely many degrees d the existence of a Pisot number of degree d and 
trace < ^i^l°^^yi/2 was proved. Theorem El improves this bound to —cy/d/ log d for some 
positive constant c. 

One ingredient needed for the proof of Theorem is a result concerning the exponent of 
maximal torsion cosets on a variety (Theorem |H)), which may be of independent interest. 
To end the introduction, we mention one immediate consequence of Theorem ^ 

Corollary 3. For infinitely many n there is a totally positive algebraic integer of degree n 
and trace less than 2n — j log log n/ log log log n. 

This follows easily from Theorem ^ using the fact that r + l/r + 2 is totally positive 
for any Salem number r. 

2. Outline of the proof 

There are two main ingredients in the proof of Theorem^ The first is a new construction 
for Salem numbers, which uses pairs of polynomials whose zeros interlace on the unit 
circle. It is an extension of the Salem number construction method used in [H], where 
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the interlacing polynomials arose from star-like trees. This new construction produces a 
polynomial of any specified negative trace that is, up to a possible cyclotomic factor, the 
minimal polynomial of a Salem number (or a reciprocal Pisot number). 

The purpose of the second ingredient is to get rid of the possibility of a cyclotomic factor, 
while at the same time bounding the degree of the Salem number. It is based on ideas of 
Schmidt |TH], and gives an upper bound for the exponent of a maximal torsion coset on a 
variety. This result is applied to a particular hypersurface to prove that the parameters in 
our polynomial construction can be chosen so that the polynomial in fact has no cyclotomic 
factor. This gives us our Salem number of the specified trace, with a bound on its degree. 



Further, the equation f{x) = 1 has real roots 71, . . . ,7^, where ai < 71 < /3i < 0:2 < 
72 < /52 ■ ■ ■ < Id-i < Pd-i <ad< Id- Also 'yd> A if and only if f{A) > 1. 

Conversely, every f{x) of the form ^ can be written in the form (Cj) for some 7 > 
and f3d-i that interlace with the aj. 

Proof. The interlacing condition for the roots easily implies (j21). Then the results follow 
immediately on applying the Intermediate Value Theorem to Ylj ^-a • 

We say that a pair of relatively prime polynomials p and q satisfy the circular interlacing 
condition if they both have real coefficients, positive leading term, and all their zeros lie on 
the unit circle, and interlace there. This last condition means that, progressing clockwise 
around the unit circle, a zero of p and a zero of q are encountered alternately. Thus p and 
q have the same degree, and neither has a multiple zero. Note too that if p and q satisfy 
the circular interlacing condition, so do p(z") and for n = 1, 2, . . . . In particular, the 
pair 2;" — 1 and + 1 satisfy it. 

Pisot numbers whose minimal polynomials are reciprocal behave in some ways like Salem 
numbers. It is clear that they must be quadratic. 

Proposition 5. Suppose that the polynomials p and q satisfy the circular interlacing con- 
dition, have integer coefficients, and that p is monic (and thus cyclotomic). Then 

(a) if piX) = 0, or q{l) = and 2p(l) — q'{l) < 0, then {z^ — l)p{z) — zq{z) is the 
minimal polynomial of a Salem number (or perhaps a reciprocal Pisot number), 
possibly multiplied by a cyclotomic polynomial. [Note: one of p(l) and q{l) is 
always 0.] 



3. Construction of Salem and Pisot numbers by interlacing 
Lemma 4. Suppose that 7 > 0, ai < /3i < 0^2 < • • ■ < Pd-i < ctd ^ A, and 




(1) 



Then f{x) can be written as 




with Xj > for all j. 



(2) 
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(b) always [z^ — z ~ ^)p{z) — zq{z) is the minimal polynomial of a Pisot number. 

Proof. Firstly, it is clear that, as the zeros of p and q interlace, both 1 and —1 must be 
zeros of pg, all other zeros of both p and q occurring in complex conjugate pairs. Put 
z + 1/ z = X, real and in [—2, 2] for z on the unit circle. 

(a) Let 7 be the leading coefficient of q. 

Suppose first that p and q have even degree 2d. If 2;^ — 1 divides g, then 

where the aj = rj + l/rj, (3j = Sj + l/sj for zeros rj of p, Sj of q, and 

-2 < ai < /3i < a2 < ■ ■ ■ < Pd-i < aa < 2. 

Thus, by Lemma IH this quotient is equal to for some Aj > 0. On the 

other hand, if 2;^ — 1 divides p, then 

z q{z) _ 711^=1(3^ -/^j) 
z'-l'p{z) - (a;2-4)n,t;(a:-a,)' 

with 

-2 < /5i < ai < /?2 < ■ ■ ■ < ad-i < (3d<2. 
Thus Lemma 0] can be applied again to give the same conclusion. 

Now suppose that p and q have odd degree 2d + 1. Then for e equal to one of 
±1, {z — e) divides q and {z + e) divides p. Then 

z q{z) _ lYY^j=i{.x - I3j) 



z^ 



1 P{z) (a; + 2£)n,ti(a;-a,)^ 



and 



-2 < I3i < ai < (32 < ■ ■ ■ < I3d < ad <2 for e = 1 

-2 < ai < I3i < a2 < ■ ■ ■ < ad < I3d < 2 for e = -1 

so that again Lemma El applies. 

Now one of {p(l),g(l)} is zero and the other positive, and (z^ — ^)p{.z) — zq{z) 
is the numerator of 1 — f{x). Hence the conditions given at z = 1 are clearly those 
necessary and sufficient for (z^ — l)p{z) — zq{z) to have a real zero greater than 1. 
(b) We consider the sum 

q{z) 2" + 1 
p{z) — 1 

and write its uncancelled numerator and denominator as q*{z)g{z) and p*{z)g{z) 
respectively, where q* and p* are relatively prime. Then, because the pairs {q,p} 
and {z"- + 1, z"' — 1} satisfy the circular interlacing condition, so do q* and p*, by 
Proposition ini (Note that there is no circularity, as the proof of Proposition [pluses 
only part (a) of this proposition.) Then g{z){{z'^ — l)p* — zq*) = z"^ {{z'^ — z — l)p — 
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zq) — {z'^ + z — l)p+zq is of the form z"- R{z) ± R* (z) , where R{z) = {z^ — z — l)p — zq 
and R*{z) = z'^'^^ ^ R{1 / z) , using the fact that one of p* and q* (say p*) is reciprocal, 
and the other (say q*) satisfies z'^'^^'^*q*{l/z) = —q*{z). Now, for any e > 0, apply 
Rouche's Theorem on the circle \z\ = 1 + e, and let n — >■ oo. This shows that 
R{z) = {z"^ — z — l)p — zq has at most one zero in 1^1 > 1. Also, R has no zeros 
on 1 2; I = 1, as any such zero would also be a zero of R*, and, as is easily checked, 
R and R* are relatively prime. Finally, because one of {p(l),g(l)} is zero and the 
other positive, in fact -R(l) < 0, so that R does have one zero on z > 1. 

□ 

Proposition 6. Suppose that the pairs of polynomials Pi^qi {i = 1,...,/) each satisfy 
the circular interlacing condition. Then is equal to a quotient where p and q 

also satisfy the circular interlacing condition. 

Proof. From Lemma 0] and the proof of Proposition El^a), we know that for each i 

z'^ — 1 Pi{z) ^ X — aj 

where x = z + 1/z, the Aj are positive, and the aj are all real and in [—2, 2]. On adding, 
the same applies to J2i ' Hence, by Lemma El this sum is equal to a positive 

scalar multiple of a quotient of polynomials t^t^, ^, where — 2 ^ ai < /3i < 0:2 < ■ ■ ■ < 

Pd-i < ttd ^ 2. Then on substituting x = z + 1/z and considering separately the cases 
when = —2 or = 2, we get the main result. 

□ 

4. The exponent of maximal torsion cosets 

As usual, let Gm denote the multiplicative group of C. An r-dimensional subtorus H of 
is a subgroup of the group = {(xi, . . . , x„) \ Xi ^ 0} where, for some r, parameters 
ti, . . . ,tr and integer matrix E = (eji)(j=i,...,r;i=i,...,n) of rank r we have Xi = t^" ■ ■ -f^". It 
is an algebraic subgroup of G^, defined by the equations {x*^ = 1 | a e A}, where the 
sl E A span the lattice of integer vectors orthogonal to the rows of E. A torsion coset is 
a translate u:H of if by a torsion point uj = {ui, . . . , Un), the Ui being roots of unity. An 
exponent of ujH is any multiple of its order as an element of the group G'^/H. A maximal 
torsion coset of a variety ^ is a torsion coset not properly contained in any other torsion 
coset in V. Results of Laurent [10, Th. 2], Bombieri and Zannier jSl Th. 2], and Schmidt 
[THl pp. 159-60] state that for any variety V C defined over a number field K, the 
union of all torsion cosets contained in V is in fact contained in a union of finitely many 
maximal torsion cosets in V, with an upper bound for this number depending only on the 
parameters of K and V. Furthermore, in jlSj Schmidt has given an explicit bound of this 
kind. 

The finiteness of the number of maximal torsion cosets in V immediately implies the 
existence of a single exponent for all these cosets. This fact can be used to prove, as in 
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Section ini that there are Salem numbers of a given trace, but without the upper bound on 
the smallest degree of such a number. The results that follow (Corollary El in particular) 
are needed to produce this degree bound. 

We denote a typical torsion coset by C = ujt^ = (i^i nj=i ^j'^^')(«=i, ■■■,«) *^ ^ being 
an r X n integer matrix of rank r. 

Consider a system of linear equations 

N 

Y^auu^ = Q (£=1,...,L). (5) 

i=l 

Following Schmidt, a solution u = (mi, . . . ,uj<i) G will be called nondegenerate if there 
is no subset / of {1, ... , A^} with < < N such that 

aaUi = (£=1,...,L). 

Lemma 7. (see [IHl p. 168-9], ^ ) Suppose we have a nondegenerate solution of ^ 
where the Ui are all roots of unity. Then, up to a factor of proportionality, the Ui are all 
P^-th roots of unity, where is the product of all primes up to N. 

In fact, their result tells us that such solutions are m-th roots of unity, where m is the 
product of at most distinct primes p ^ N. However, we need an exponent valid 

uniformly for solution sets of different such A^-term equations. This is why we take -th 
roots of unity, being the 1cm of all such m. A uniform 'killer' exponent is provided by 
the following result, and its corollary. 

Theorem 8. Suppose that V is an affine variety in GJ^ defined over Q, given say by 
polynomial equations 

^a«x' = (£=1,...,L) (6) 

i 

with total degree d. Suppose also that the set 

MiV) = {i G Z" I a« ^ for some i} 

has diameter D{V). Then every {n — k)- dimensional maximal torsion coset on V has an 
exponent mPj^ for some integer D{Vf^k''l'^. Here N := #7V(V) ^ ("+^) . 

Proof. The ingredients for the proof come from Schmidt ^Hl- Take r = n—k and a maximal 
r-dimensional torsion coset C = ujt^ on y, so that 

5^a,ia;'t^^ = (£=1,...,L). 

Our aim is to find uji with also C = uJit^, with cji a vector of {mP^ )-th roots of unity for 
some m ^ DiyY^k^^'^ . Now for any j G Z'' the coefficient of is 



5^a«a;^ = (£=1,...,L). 

i:Si=j 



(7) 
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Here the sums over i are taken over all relevant i in MiV). Now (j?)) may be degenerate, 
splitting into nondegenerate equations 

5^a,ia;' = (£ = 1, . . . , L, g G Q say) (8) 

for nonempty subsets Iq of Z". Now, for a single g, apply Lemma [3 to (jHI), to obtain 

5^a,ia;'-'« = (^=1,...,L) (9) 

where \q is some fixed vector in Iq. Here, the number of terms is at most . Then for all 
q E Q, we have from Lemma d that all uj^~^''{q e Q) are vectors of i^-th roots of unity. 

Recalling that k = n — r, we claim that the set of all vectors {i — ig | i G Iq,q G Q} 
generates a fc-dimensional sublattice Cc of Z". For, from ((7j), the lattice spanned by the 
rows of E is orthogonal to Cc, and so Cc has dimension ^ n — r. But if the inequality were 
strict, there would be a nonzero vector i' G Z" orthogonal to Cc and not in the rational 
span of C^. Then for i G we would have i' ■ i = i' ■ ig (and also of course Ei = Eiq), so 
for any u G we would have, for i = 1, . . . , L, 

<? ieiq 
= 0, 

and so the larger torsion coset cut'^ would lie on V, where t' = (t,^) and E' = (^), 
contradicting the maximality of ut^. 

Next take a basis of vectors in {i — ig | i G Jg, g G Q} for Cc, and put 

u; = e*^ = (e*^i, . . . , e'^"). Write = Y!1=i ^j'^j + V', where t}) = pE for some p 
Then, on solving the system of linear equations 

k 

£i-e = Y^ X/i -Ij {i = l,...,k) 

and using the fact that Pj^ti ■ 6 = (mod 27r) {i = 1, . . . ,k), we see that det(£j ■ 
ij)X = (mod 27r). Note too that det{ii ■ £j) ^ 0. Then, using the Cauchy-Schwartz and 
Hadamard inequahties, we have that 

|det(£i-£,)| ^ Diyf^k''!''. 

Put ti = e"*''. Then tf = e"*"^, and for lc\ = ujtf and some m' = mP^ with m ^ 
D{Vf^k''l'^ we have a;™' = 1. Since 

C = Ljt^ = Lj(tit)^ = LJita^ 

say, we see that C has exponent m' . □ 
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This result immediately gives us a killer exponent K valid for all maximal torsion cosets 
on V . 

Corollary 9. Let V he as in the Theorem, and K = Pj^ lcm(l, 2, . . . , D{V)'^^n'^^'^), where 
N = ^Miy). Then every maximal torsion coset ofV has exponent K. 

5. Maximal torsion cosets on a particular hypersurface 

We shall be applying the results of the previous section to the affine hypersurface h{x) = 

0, where x = {xq, . . . , Xn) G and 

n n n 

/i(x) = 2{xl - 1) llixi - 1) - Xo ^(x, + 1) l[{xi - 1). 

i=l j=l i=l 

The reason for looking at this hypersurface is that we shall apply the identity 

^(x) _ 3^0 - 1 Is^ Xj + l 

2a;onr=i(3^^-l) ^0 2^x,-l' ^ ' 

which is connected to our interlacing considerations of Section El 

Lemma 10. The only maximal torsion cosets of h with xq nonconstant are the algebraic 
subgroups Bij of Q"^^ , where i ^ j are both nonzero, and 

Bij = {x I Xi = Xj = l,xo = to,xe = te{i 7^ hj)}, 

of rank n — 1. 

Proof. Clearly no point on h = can have just one Xj = 1. If x with any two Xj,Xj = 
1 is on /i = then it belongs to Bij. Thus any other rank r maximal torsion coset 
with Xq nonconstant has no Xj identically 1, so that we must have Xq = ujQtl^° ■ ■ ■ t'^''° , 
where (eio, . . . , e^o) 7^ and Xj = ujitl^' ■ ■ ■ t'p'\ where (eij, . . . , e^) 7^ whenever oji = 

1. By avoiding certain hyperplanes we can choose ±(fci, . . . ,kr) G Z'*" not orthogonal to 
any nonzero (ch, . . . , Then for {ti, . . . ,tr) = {t^^ , . . . ,t^''), Xi = Uit^'- where ii := 
Yl]=i ^j^ji 7^ when Ui = 1, and, by choice of the sign, io > 0. Now as t 00, the 
right-hand side of (jlUj) goes to infinity, so that the coset cannot be on /i = 0. □ 

We now estimate the killer exponent K, valid for every maximal torsion coset on this 
hypersurface, defined over G^^. 

Lemma 11. There is a killer exponent K with log log < 0.2 + (3(n + l)/2) log(n + 3) 
for the hypersurface h = 0. Further, K can be chosen with all its prime factors less than 
(n + 3)3("+i)/2. 

Proof. The hypersurface has diameter D = ^/n + 4, degree d = n + 2, N = ^Af{h = 
0) = 3 ■ 2". Hence D^'^+^in + l)("+i)/2 < (n + 3)3(n+i)/2^ Corollary H gives K = 
Pj^ ■ lcm(l, 2, . . . , D^''+\n + l)(«+i)/2)^ with all prime factors of K less than {n + 3)3(n+i)/2, 
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Then, using standard bounds of Rosser and Schoenfeld JHl for the arithmetical functions 
djip we obtain 

logK < 9{3 ■ 2") + + S)^^"^^)/^) 

< 1.02 • 3 • 2" + 1.04 ■ (n + 3)=^("+^)/2 

< 1.2(n + 3)=^("+l)/^ 

giving the upper bound claimed. □ 

6. Proof of Theorem [T] 

The following lemma will complete the proof of Theorem ^ 

Lemma 12. For given even n there are positive integers ki, . . . ,kn such that 

/i(t,t^\...,t'=") = 2(t- l)"-i5(t), 

where S{t) G is monic irreducible and the minimal polynomial of a Salem number of 
trace T := 1 — n/2. Further, S has degree less than expexp(22 + 4TlogT). 

Proof. For all maximal torsion cosets of h with xq constant (ie all except the Bij) we can 
suppose that the constant Xo-values are all K-th roots of unity, where furthermore K has 
been chosen minimally. Note that K is certainly even, because the point Xq — Xi — • • • — 
Xj\f = —1 lies on h = and, as it is on no B^j, must lie in one of the constant-Xo maximal 
torsion cosets. Take ki = K, and k2, . . . ,kn as the smallest n — 1 primes not dividing 
K. Then all fci, . . . , fc„ are pairwise relatively prime. We now assert that for every root of 
unity uj and uj^ = {uj, uj^^, . . . , cj^") with h{uj^) = 0, we have u = 1. 

For uj^ belongs to some maximal torsion coset. If uj^ has at least two components = 1, 
then (by the extended euclidean algorithm) uj = 1. Alternatively, it belongs to no Bij, and 
so to some maximal torsion coset with xq constant, Xq = u, and so = uj^^ = 1. This is 
impossible, as we cannot have just one Xj = 1, as noted above. This proves the assertion. 

It is easy to check that {d/dt)"'h{t,t''^ , . . . ,t^") evaluated at t = 1, is nonzero. Further- 
more, h{t, t^\...,t^^) = nr=i(^i - l)(2(t^ - 1) - tn) = (mod 2), so that all coefficients 
of h are even. This gives the stated factorization 2(t — l)"^^S'(t) of h{t, t'^i, . . . , t'^"). Also, 
3jS clII ki ^ 2 , 

h{t,t''\...,t^") = 2t2+Sfc. + 

showing that S{t) has trace 1 — n/2. 

Finally, to show that 5* is the minimal polynomial of a Salem number, observe that we 
have shown that none of its zeros are roots of unity. Now since t^ + 1 and — 1 satisfy the 
circular interlacing condition, so does the sum | ^"^^ Proposition IHl so we can 

write it as q{t)/p{t), where p and q satisfy the circular interlacing condition. Furthermore, 
as n is even, q has integer coefficients, as does p = Ylii^''' ~ ~ 1)""^. Hence, as 
p{l) = 0, the numerator S{t) = {t^ — l)p{t) — tq{t) of the right-hand side of (fTUj) with 
Xq = t,Xi = t^^ is, by Proposition Efa), the minimal polynomial of a Salem number of trace 
1 - n/2. 
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Now the degree of 5* is 2 + J2i h — (n — 1), and from Lemma [TT] we can take k2, . . . ,kn 
to be the smallest n — 1 primes greater than {n + S)^*^""'"^-*/^. By Bertrand's Postulate 
(Chebyshev's Theorem), this gives degS < K +{n- l)2("-i) ■ {n + 3)3("+i)/2 < 2K, and 
loglogdegS < loglogir + log2/log/s: < 0.2 + (3(n+ l)/2) log(n + 3) +0.1. For n = 2T + 2 
one readily checks that this is less than 22 + 4T log T. □ 

Remark 13. There are many maximal torsion cosets with Xq constant, for instance Xq = 
—l,Xi = a;^^ = ti, . . . ,Xn-i = = tn/2- Also one can for instance construct some for 
xo = 1 using the identity 3 cot n/S — cot 7r/6 = 0. 



7. Proof of Theorem [51 

Proof. The proof is much easier for Pisot numbers, as there are no possible cyclotomic 
factors to dispose of. We replace the fraction (t^ — l)/t in (fTTH) by (t^ — t — l)/t, and 
then can simply choose the parameters ki to be the first n primes. Thus, again putting 
\ T^l=i = Q{'t)/p{t), the polynomial {z'^ — z — l)p{z) — zq{z) will be the minimal 
polynomial of a Pisot number of trace 2 — n/2. □ 



8. Computing Salem and Pisot numbers of negative trace. 

Salem and Pisot numbers of negative trace can be produced using as m 

the previous sections. Thus, for the Pisot numbers of trace — T, the first 2T + 4 primes are 
used for the k^. For the Salem numbers of trace — T, the first 2T + 2 primes are used for 
the ki. In particular, ki is taken to be simply 2, instead of the very large killer exponent K 
used in the proof above. Computation using Maple shows that this produces a polynomial 
free of cyclotomic factors for T ^ 25, giving a Salem number of trace — T and degree equal 
to the sum of the first 2T + 2 primes minus 2T — 1 (for instance degree 5540 for trace —25). 
However, we do not know whether this always happens. It would of course be nice if this 
could be proved, as we would then obtain a degree bound in Theorem [T] as good as that 
in Theorem |2l 

Here is some pseudocode that gives the minimal polynomials. For a Salem number of 
trace — T: 



r = l;S={z^-l){z-iy,Q = z; 

for j = 1, . . . , T + 1 do 

q = nextprime(r); r = nextprime(g) ; 



enddo 

if gcd{S{z),S{-z)S{z^)S{-z^)) = 1 then print(5'); 
endif 



The gcd condition tests for cyclotomic factors, and is based on the fact that a root of 
unity u is conjugate to one of —u, or —u^. See for further developments of this idea. 
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For instance, for T = 2 we obtain the (reciprocal) Salem polynomial 

S-2{z) = z^^ + 2 z^'^ - 2 - 19 z^^ - 57 z^* - 123 z'^^ - 222 z^^ _ 357 ^31 _ 537 ^,30 _ 727 z'^^ - 950 2^** - 1190 z- 



Note that the Salem polynomials S produced by this method have |S'(— 1)S'(1)| large. 
This is easily seen by putting z = ±1 in the pseudocode. An interesting question is 
whether there are Salem numbers with arbitrary trace and |S'(— 1)S'(1)| = 1, the so-called 
unramified Salem numbers (see Gross and McMuUen |8j). 

For a Pisot number of trace — T: 

r = l-p = z"^ - z-l-Q = z- 

for j = 1, . . . ,T + 2 do 

q = nextprime(r); r = nextprime(g) ; 

P ={z'i - 1)(;2^ -I). P- (^e+r- - 1) . g; 

enddo 
print(P); 

Finally, we justify the statements in the Introduction. We note that there are Salem 
numbers of every nonnegative trace: for > the polynomial z^—nz^ — {2n+l)z'^—nz+l = 
z'^{{z + l/z)'^ — n{z + l/z) — (2n + 3)) is easily seen to be the minimal polynomial of a Salem 
number of trace n. This follows from the fact that x"^ — nx — {2n + 3) has one zero in (—2,2) 
and the other zero greater than 2. Also, z^ — z"^ — 2z^ — z"^ = z^{{z + l/z)^—A{z+l/z)—2) 
is the minimal polynomial of a Salem number of zero trace. 

The lower bound [1 — d/9\ for the trace of a degree d ^ 10 Salem number follows from 
the fact that the trace of a totally positive algebraic integer of degree n ^ 5 is greater than 
16n/9 (^21) on noting that for a Salem number r of degree d and trace — T, the number 
r + l/r + 2is totally positive of degree d/2 and trace d — T. Thus (turning the inequality 
around) for — T ^ —2 every Salem number of trace — T has degree at least 



This inequality is sharp for T = 2 (|I31)- The only Salem number of degree less than 10 
having negative trace is the one with minimal polynomial z^ + z'^ — z^ — 4z^ — 5z'^ — 4z^ — 
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